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for the order. 

As a consequence of our approach we deduce an upper bound for the Hilbert- 
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implicit differential systems. 

Key words: Implicit systems of differential equations, differentiation index, order of a 
differential ideal, existence and uniqueness, differential membership problem. 

*Partially supported by UBACyT X112 (2004-2007) and X211 (2008-2010). 

tPartially supported by CONICET PIP 5852/05, UBACyT X847 (2006-2009) and ANPCyT PICT 
17-33018/05. 

IPartially supported by UBACyT X058, CONICET PIP 5056, ANPCyT PICT 03-15033 and Fundacion 
Ciencias Exactas y Naturales (Argentina). 



1 



1 Introduction 



Throughout the paper we consider imphcit differential algebraic equation systems (DAE 
systems for short) of the following type: 

Jl\^l, ^1, ■ ■ ■ , , ■ ■ ■ , ^n, ^n, ■ ■ ■ , ) — U 

(s) := ; (1) 

f t Y Y yC^i-l) Y Y v(ern)\ r, 

V Jr\-^1, ^1, ■ ■ ■ , , ■ ■ ■ , y^n, ^n, ■ ■ ■ , ) — ^ 

where r < n, and /i, . . . , are polynomials in the n differential unknowns X := Xi, . . . , Xn 
and some of their derivatives with coefficients in a characteristic zero differential field K 
(i.e. a field with a derivation, for instance := Q,M or C with the null derivation, or a 
field of rational functions in one variable with the standard derivation). Each non- negative 
integer e^j, 1 < i < r, 1 < j < n, denotes the maximal order of derivation of the variable 
Xj appearing in the polynomial /j. 

The set of all total formal derivatives of the equations defining (S) generates an ideal 
[F], closed by differentiation, in the (algebraically non-noetherian) polynomial ring K{X} 
in the infinitely many variables j = 1, . . . ,n, I £ Nq. Since this ideal is not necessarily 
prime, we also fix a minimal prime ideal ^ C K{X} containing [F]. 

We require a regularity assumption on the system (E) with respect to the prime ideal 

we suppose that the system itself, as well as the systems obtained by successive to- 
tal derivations of the defining equations, are non-singular complete intersection algebraic 
varieties at the (not necessarily closed) point ^p. In terms of Kahler differentials this con- 
dition is formalized saying that the differentials {dfj^^\ I < i < r, k G No} C ^k{x}\k 
are a linearly independent set in ^k{x}\k ^ DAE system 

with this property is called quasi-regular at *P (see, for instance, |21j . [31j). In addition, 
we make use of a rather technical hypothesis which holds for a wide class of DAE systems 
(see Subsection 12.41 below) 

The main invariant we consider in this paper is the differentiation index of the system 
(S). There are several definitions of differentiation indices not all completely equivalent 
(see for instance [I], [3], [11], [M], [29], [33], [M], [36], [37], [38], [53]), but in every case 
it represents a measure of the implicitness of the given system. For instance, for first 
order systems, differentiation indices provide bounds for the number of total derivatives 
of the system needed in order to obtain an equivalent ODE system (see [U Definition 
2.2.2]). Thus, differentiation indices turn out to be closely related to the complexity of 
the traditional numerical methods applied to solve these systems (see [H Theorem 5.4.1]). 
Since explicitness is strongly related to the existence of classical solutions, a differentiation 
index should also bound the number of derivatives needed in order to obtain existence and 
uniqueness theorems (see [35], [36], [37]). 

Following the construction given in our previous paper [H], we introduce here a new 
differentiation index for a DAE system (S) quasi-regular at a minimal prime ideal ^p, 
that we call the "^-differentiation index. As usual, its definition follows from a certain 
chain which eventually becomes stationary. In our case, this chain condition is simply 
established by the sequence of ranks of certain Jacobian submatrices associated to the 
input equations and their total derivatives (see Theorem [8] below). This approach enables 
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us to show in an easy way several properties and consequences of our differentiation index. 

In particular, we show that the *P-differentiation index is closely related to the number 
of derivatives of the system needed to obtain the manifold of all constraints that must be 
satisfied by its solutions. More precisely, we prove that for every order h, all the differential 
consequences of order h of the quasi-regular system (S) can be obtained from the first 
h + a — maxjejj} derivatives of the equations (see Theorem [TUj) . 

The second invariant of the system (S) we consider is the order of the prime differential 
ideal ^ (which will also be called the order of (S) in the case where the differential ideal 
generated by the equations is prime). Roughly speaking, the order of a DAE system is 
the number of initial conditions that can be prefixed arbitrarily. 

In his posthumous papers [16] and [TT], Jacobi introduced the parameter J{£) := 
max{^[^^ ej^(j) | r : {1, . . . , r} ^ {1, . . . , n} is an injection} as an optimal estimation for 
the order of an ordinary DAE system (here eij := —oo if the variable Xj does not appear 
in fi). This result should be considered as a conjecture, since Jacobi's proof seems not to 
be complete in the general case. Nevertheless, in the zero-dimensional case, it has been 
proved for linear systems by Ritt [lOl Ch. VII, p. 135] and extended to quasi-regular DAE 
systems in (see also [S])- Other less accurate upper bounds for the order, mainly of 
Bezout type, have been given for general (zero-dimensional) systems (see, for instance, [6], 

[H], m Ch. IV, Prop. 9], m, m)- 

Here, we establish an a priori upper bound for the sum of the ^^-differentiation index 
and the order of in the quasi-regular setting, in terms of the maximal derivation orders €ij 
of the variables in the input equations. More precisely, we prove the following inequality: 

a + ord(*p) < J (So) + max{eij} - min{ei-,}, 

where a denotes the *P-differentiation index of the system (S), ord(*p) the order of the 
differential ideal *P and J(fo) the weak Jacobi number J(fo) := max{^[^j^ ^iT{i) I '■ 
{1, . . . , r} —5- {1, . . . , n} is an injection} with €ij := if the variable Xj does not appear in 
fi. In particular, we deduce an upper bound for the differentiation index in terms of the 
weak Jacobi number J{£o) (see [8], [35], [37] for less precise estimations). 

Moreover, our techniques enable us to recover the original Jacobi bound for the order 
of quasi-regular DAE systems (see Theorem llSp . Our strategy may be considered as a 
generalization of Lando's method introduced to obtain a weak Jacobi bound for the order 
of linear differential equations (see [28] ) . 

Our approach, which is reminiscent of the classical completion techniques in partial 
differential algebraic equations (see [H], [IB], [27] and [33), also enables us to exhibit 
an upper bound for the Hilbert-Kolchin regularity of the ideal ^ depending only on the 
maximal order of derivation of the variables involved in the system (see Theorem I12|) . A 
precise estimation of this regularity can be obtained if a characteristic set of the ideal ^ 
is known (see \20\ Ch.II, Section 12, Th.6 (d)]) while in our case this requirement is not 
necessary. 

As a further consequence of the previous results, we deduce an effective differential ideal 
membership test for quasi-regular DAE systems (Theorem I20p : we obtain asymptotically 
optimal bounds on the order and degrees for the representation of a differential polynomial 
lying in the ideal in terms of the given generators (see [12], [2] and [13] for related works). 
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Finally, in the zero-dimensional case (n = r), we obtain a result concerning the number 
of derivatives of the input equations required to obtain an explicit ODE system (see 
Theorem \22\i and an existence and uniqueness theorem (see Theorem I24p . 

The paper is organized as follows: in Section [5] the notion of a quasi-regular DAE 
polynomial system is introduced, along with some basic notions from differential algebra 
we use. In Section [3] we give the precise definition of the *P-differentiation index and 
show the relationship of this invariant with the manifold of constrains of the system. The 
Hilbert-Kolchin regularity and the order of the ideal ^ are analyzed in Section [H Section 
El is devoted to the proofs of Jacobi-type bounds for the differentiation index and the order. 
We present our results on the differential membership problem for quasi-regular systems 
in Section [6l The last section of the paper is concerned with existence and uniqueness of 
solution results for implicit quasi-regular DAE systems. 



2 Preliminaries 

2.1 Basic Notations 

Let K he a characteristic zero field equipped with a derivation 5. For instance K = 
or C with 6 := 0, OT K = Q{t) with the usual derivation 5{t) = 1, etc. 

For an arbitrary set of (differential) indeterminates Zi, . . . , Za over K we denote the 
p-th successive derivative of a variable Zj as zj^^ (as customarily, the first derivatives are 

also denoted by Zj); we write Z^f) := {z[^\ . . . , Z^^} and Z^ := {Z(^), 0<i<p}. 

The (non-noetherian) polynomial ring K[Z^p\ p € Nq], called the ring of differential 
polynomials, is denoted by K{Zi, . . . , Zq,} (or simply K{Z}). 

Given a finite set of (differential) polynomials Hi, . . . , Hp G K{Zi, . . . , Za}, we write 
[Hi, . . . , Hp] to denote the smallest ideal of K{Zi, . . . , Za} stable under differentiation, 
i.e. the smallest ideal containing Hi, . . . ,Hp and all their derivatives of arbitrary order. 
The ideal [Hi, . . . , i?^] is called the differential ideal generated by Hi, . . . , Hp. 

For any differential polynomial g lying in a differential polynomial ring K{Zi, . . . , Za} 
the following recursive relations hold for the successive total derivatives of g: 

9^'^ := 9, 

where 6{g^~^^) denotes the polynomial obtained from g^P~^^ by applying the derivative 6 
to all its coefficients (for instance, if K is a field of constants, this term is always zero). 

2.2 The system 

Let r < n € N. Throughout the paper we consider DAE systems of the following type: 

f /i(xf"i,...,xH) = 



(2) 



rl^l 5 • • • ! 
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where fi, . . . , fr are differential polynomials in the n differential variables X := Xi, . . . , Xn 
with coefficients in the field K. Each non- negative integer €ij denotes the maximal deriva- 
tion order of the variable Xj appearing in the polynomial fi. We denote e := maxjeij} for 
the maximal derivation order which occurs in (S). We assume that (S) actually involves 
derivatives, i.e. e > 1. 

We write [F] C K{X} for the differential ideal generated by the polynomials fi, . . . , fr. 

We introduce also the following auxiliary (noetherian) polynomial rings and ideals: for 
every i S Nq, Ai denotes the polynomial ring Ai := and Aj C ^j_i+e the ideal 

generated by the total derivatives of the defining equations up to order i — 1, namely 
Aj := (/|* ^\ . . . , fr (this ideal is usually known as the i — 1 prolongation ideal). We 
set Aq := (0) by definition. 

2.3 Quasi-regular DAE systems 

The notion of quasi-regularity appears implicitly in p9] in order to generalize a Janet 
Conjecture to non-linear systems. Up to our knowledge it is the more general frame where 
the Jacobi order bound holds (see for instance [7], [21], |31j . |32j). 

Definition 1 C K{X} is a prime differential ideal containing fi, ■ ■ ■ , fr (or equiva- 
lently, containing [F]) we say that the system (S) is quasi-regular at *p if for every positive 

integer i, the Jacobian matrix of the polynomials /|* ^\ . . . , fr with respect to the set 
of variables has full row rank over the domain ^j_i4.e/(^j_i+e H ^). 

This condition can be easily rephrased in terms of Kahler differentials (as in Johnson's 
original work [19]) saying that the differentials {d/j*\ 1 < i < i ^ Nq} C ^k{x}/k 
a linearly independent set in ^k{x}/k ^k{x} 

Geometrically, Definition [1] means that for any positive integer i the algebraic variety 
Vi given by the ideal generated by the ri polynomials /j* . . . , fr in the (i — 1 + e)n- 
variate polynomial ring Aj„i_|_e is smooth in almost all point of the closed subvariety 
defined by the prime ideal Ai_i+e n *p. 

This notion can also be interpreted in the context of jet spaces: for the sake of sim- 
plicity, assume that ^ is the prime ideal associated with a solution of the ideal [F] having 
infinitely many coordinates. The quasi-regularity condition means that, if ip{t) is a so- 
lution of the differential system such that (p^^\to) € M" are the j-th. coordinates of the 
point defined by ^p, for each integer i the jet of (p at tQ is a regular point of the subspace 
V C J°°(M,M") defined by the input system differentiated up to order i. 

In the sequel, we will assume that not only a system (S) is given, but also a prime 
differential ideal *p C K{X} minimal with the property of containing the polynomials 
fi, . . . , fr and such that (S) is quasi-regular at *p. 

For each non-negative integer i we write Bi for the local ring obtained from Ai after 
localization at the prime ideal AiD^. Since each Ai is a polynomial ring, the localizations 
Bi are regular rings (see for instance \26\ §1, Corollary 1.8]). We denote ^i the maximal 
ideal generated by n ^ in B^ . For the sake of simplicity we preserve the notation Aj 
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for the ideal generated by the derivatives up to order i — 1 of the polynomials /i , . . . , 
in the ring Bi^i+e- 

Remark 2 If the ideal [F] C K{X} is already a prime ideal, the minimality o/*p implies 
that *p = [F] and all our results remain true considering the rings Ai and the ideal [F] 
without localization. In this case if (S) is quasi-regular at [F] we will say simply that (S) 
is quasi-regular. 

We mention now some easy consequences of the quasi-regularity hypothesis related to 
localization at the prime 

Proposition 3 Let (S) he a DAE system defined by polynomials F := fi, . . . , fr whose 
orders are bounded by e and let *P C K{X} be a minimal prime differential ideal containing 
F such that (S) is quasi-regular at Let i (^N be an arbitrary positive integer. Then we 
have: 

1. /|* . . . , fl- '^^ is a regular sequence in the local ring i?j_i+e CLnd generates a prime 
ideal. 

2. In the localized ring K{X}s:^, the polynomials fi,. . . ,fr generate the maximal ideal 
VK{X}<^. 

3. // IK denotes the residual field of the prime ideal the differential transcendence 
degree of IK over K is n — r. 

Proof. The first assertion is a direct consequence of the classical Jacobian Criterion 
(see for instance [301 §29]) applied to the regular local ring i?j_i+e and the polynomi- 
als /j* ^\ ■ ■ ■ , fr which belong to the maximal ideal ^j_i+e. So, the quasi-regularity 
implies that these polynomials form part of a system of parameters in i?j_i+e and in 
particular they generate a smooth (hence prime) ideal and form a regular sequence. 

Since only a finite number of variables are involved in primality checking in K{X}<;p, 
the previous arguments imply that the differential ideal [F] K{X}(;p is necessarily prime 
in the local ring K{X}fp and then, by minimality, it agrees with *pi^{X}(p. 

For the last assertion we observe first that the polynomials /i , • • • , /r are in fact dif- 
ferentially algebraically independent over K: any finite family of them and their deriva- 
tives is part of a system of parameters of a suitable local X-algebra of type Bi for all 
i big enough, then algebraically independent over K. Now, let 6i, . . . ,9n-r G ^{^} 
be such that fi, . . . , fr,Oi, . . . , 6n-r £ K{^} is a differential transcendence basis of the 
fraction field of K{X} over K. Taking the classes of ^i, . . . ,9n-r hi the fraction field of 
(i('{X}/[F])<p one obtains a differentially algebraically independent family. Hence, the 
third assertion follows from the second condition of the Proposition. ■ 

Prom now on we will write A for the ideal [F] K{X}<;p = ^ K{X}f;p. 
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2.4 Pseudo-Jacobian matrices 



We introduce a family of pseudo-Jacobian matrices we need in order to introduce the 
notion of differentiation index in the next section: 



Definition 4 For each /c G N and i S N>e-i (i-e. i G Z and i > e — 1), we define the 
kr X kn-matrix Zk i o,s follows: 



-e+l) 







Qp(i-e+2) Qp(i-e+2) 






~dX>J+W 



\ 



where each g^fq) denotes the r x n— block consisting in the Jacobian matrix of the polyno- 



mials f[^\ . . . , fr^^ with respect to the variables x['^\ 



(9) 



Observe that the block triangular form of Zk,i follows from the fact that the differential 
polynomials have order bounded hy i + p. Hence, their derivatives with respect 

to the variables X^*"'"-') are identically zero for j > p + 1. 

Definition H] is quite ambiguous because this matrix may be considered over different 
rings, and so, invariants as the rank or the solution space are not well defined and in fact 
they may differ. 

We add the following last hypothesis on our input system (E): 



Hypothesis We assume that for any pair of indices i, k the rank of the matrix 2k,i 
over the integral domain -Bj+fc+s/^i-i+e+fc+s does not depend on s. In other words, 
the rank of each matrix Zk,i considered alternatively over the rings Bij^k/ ^i-i+e+k^ or 
Ai^k/{^i+k n ^), or the residual field of A+fc, or the field IC := K{X}/^i, is always the 
same. 



This hypothesis is satisfied for relevant classes of DAE's: for instance, linear or linear 
time varying DAE's, or DAE systems with generic second member coming from Control 
Theory (see [H] ) . This kind of assumptions can be regarded as an algebraic counterpart of 
usual hypotheses concerning solvability, smoothness and constant rank properties for the 
system and its prolongations, required in various notions of differentiation indices in the 
analytic-numerical framework (see for instance [3] or [1]). 



3 A notion of ^—differentiation index 

3.1 A linear algebra— based definition 

In this section we introduce a notion of a differentiation index of the system (S) depend- 
ing on the chosen minimal prime differential ideal ^ containing it, always assuming that 
the system is quasi-regular at this prime ideal and verifies the hypothesis introduced in 
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Subsection [231 

We introduce a double sequence of non-negative integers associated with the ma- 
trices ^k,i- 

Definition 5 For G No and i G N>e_i, we define ^k,i £ No as follows: 

- fiQ^i := 0. 

- IJ'k,i '■= diniK ker(5^ J, for k > 1, where 5^ ■ denotes the usual transpose of the matrix 

In particular Hk,i = kr - rankK(3fc,j). 

The sequence fik,i is strongly related with some algebraic facts concerning the algebraic 
dimension of the ideals Ap generated by the first (p — 1)— th total derivatives of the 
polynomials /i, . . . ,/,.: 

Proposition 6 Let A; G No and i G N>e-i. Then: 
(i) The transcendence degree of the field extension 

Prac(Si/(Ai_e+i+fc n Bi)) ^ Fmc{Bi+k/ ^i-e+i+k) 
is k{n -r) + fik,i- 
(a) The following identity holds: 

trdeg^Frac(5i/(Aj_e+i+fe n Bi)) = (n - r){i + 1) + er - fik,i. 

Proof. The proof follows [HI Proposition 2 and Remark 3]. 
Consider the following diagram of field extensions: 

Frac(5i/(Ai_e+i+fe n Bi)) ^ Frac(Sj+fc/Aj_e+i+fc) 

K 

We can consider the fraction field of the ring Bi^k/^i-e+i+k as the fraction field 
of n[X'-^+^\. . .,X^'+^^]/{F^'\. . . where 9. is the fraction field of the integral 

domain i?i/(Ai_e+i+fc H Bi). By the classical Jacobian Criterion (see for instance [301 
Chapter 10, §27]), this transcendence degree equals kn — rank Zk,i = kn — {kr — ^k,i) = 
k{n -r)+ ij,k,i. 

In order to prove the second assertion it suffices to compute the transcendence degree 
of the fraction field of Bi^k/ \~e+i+k over the ground field K. 

Set p := n *P and let q C p C ^j+fc be the prime ideal such that Aj_e+i+fc-Bj+fc = 
qi?j_l_fc (from the first assertion of Proposition [3] this prime ideal q exists and it is unique 
verifying these requirements). Therefore, the transcendence degree of the field extension 
K ^ Frac(i?j_|_fc/Aj_e+i4-fc) is equal to the Krull dimension of the domain ^j+^/q and so, 
by the catenarity of the polynomial ring Ai^k, equals dim Aj+fc — dim(^j_|_fc)q. Again by the 
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catenarity and the fact that q C p we conclude dim(Aj+fc)q = dim(Aj+fc)p — dim(Aj_|_fc/q)p. 
Since /j* '^+'^1^ . . . , fr ^^^^ form a regular sequence in -Bj+fc (Proposition [3|) we infer that 
dim(Ai^k)p — dim(Ai^k/(\)p agrees with the number of such polynomials, i.e. with r(i — 
e + 1 + k). Hence we obtain: 

trdeg;^Frac(Si+fc/Aj_e+i+fc) = dim Aj+fc - r{i - e + I + k) = {i + k + l){n - r) + er. 

The results follows directly from the additiveness of the transcendence degree. ■ 

By means of linear algebra arguments one can show as in [HI Section 3] some other 
relevant properties of the sequence /Xfc j; in particular we have the following results which 
allow us to introduce the notion of differentiation index of (S) related the quasi-regular 
prime differential ideal 

Proposition 7 Proposition 11]) The sequence does not depend on the index i. ■ 
In the sequel we write simply instead of ^k,i- 

Theorem 8 ([U Theorem 9 and Corollary 14]) There exists u G Nq such that Hk < /^fc+i 
for all k < a and = l^k+i for all k > a. m 

Definition 9 The constant o" E No appearing in Theorem\^will he called the *P-differentiation 
index of the DAE system (S). // the differential ideal generated by the equations of (S) is 
itself a prime ideal, we say simply the differentiation index of (S). 

While the dependence of a on the prime ideal *p is obvious from the construction, the 
name "differentiation index" would seem rather obscure at this point. This assignation will 
be justified throughout the paper by proving that a verifies some properties customarily 
associated to the differentiation index such as estimating the number of derivatives needed 
either to obtain the set of constraints that every solution of the system must satisfy 
(Subsection 13. 2p or to convert the DAE system (S) into an explicit ODE (Subsection 17. ip . 

3.2 Manifold of constraints 

A remarkable property associated with most differentiation indices is that they provide 
an upper bound for the number of derivatives of the system needed to obtain all the 
constraints that must be satisfied by the solutions of the system (see for instance |11|). 

For instance, suppose that the input system is of first order and we look for the purely 
polynomial conditions in . . . , X„] (i.e. no derivatives of the variables appear) which 

must be verified by any solution. These conditions define an algebraic variety which is 
usually called the manifold of constraints of the system. Clearly this variety is defined by 
the ideal [F] f] K[Xi, . . . , Xn]. From the noetherianity of the polynomial ring, there exists 
the smallest non negative integer s such that the elements of this ideal may be written 
using at most s derivatives of the input equations. Under our conditions of Section O this 
minimum integer is exactly the differentiation index: 
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Theorem 10 Let a £ 'Nq be the ^-differentiation index of the system Then, for 

every i G N>e-i, the equality of ideals 



^i-e+i+a nB, = An B, 

holds in the ring Bi. Furthermore, for every i € N>e-i, the "^-differentiation index a 
verifies: a = mm{h G Nq : Ai^e+i+h ri Bi = Ad Bi}. 

Proof. Fix an index i G N>e-i. Let us consider the increasing chain (Ai_e+i+fc ni?i)fc6No 
of prime ideals in the ring Bi. From Proposition El for every k G Nq, we have 



Since fik is stationary for k > a (Theorem[H|), all the prime ideals Aj_e+i+fc ni?j have the 
same dimension for k > a and the chain of prime ideals becomes stationary for k > a. 

It only remains to prove that the largest ideal of the chain coincides with ACiBi. One 
inclusion is obvious. For the other, let / be an arbitrary element of A n Bi, then there 
exist differential polynomials h,aij G k{X}, h such that 



If N is an integer greater than i + a and than the maximal order of the variables X 
appearing in this equality, we have f £ Bi and / G Atv-c+i C B^. Then / G A^^^+i^Bi 
and, since the above chain of ideals is stationary, / G Aj_e+i+o- H B^. This finishes the 
proof of the first assertion of the Theorem. 

In order to prove the second part of the statement, for each i G N>e_i, let hi be the 
smallest non- negative integer such that Aj_e+i+/i. fi i?j = A n -B,,. By the definition of /ij, 
the transcendence degrees trdeg^(Frac(i?i/Aj_e+i+fc fl Bi)) coincide for k > hi, and so, 
fik is constant for k > hi (see identity ([3]) above). This implies that a < hi. The equality 
follows from the first part of the statement and the minimality oi hi. ■ 

Remark 11 Taking i = e — 1 in the last assertion of Theorem \10[ we have the following 
alternative definition of the ^-differentiation index: 

a = mm{h G Nq : A^, n Be-i = AD B^-i}. 

4 The Hilbert-Kolchin polynomial of the ideal ^ 

Let ^ C K{X} be a minimal prime differential ideal such that the system (S) is quasi- 
regular at The last assertion in Proposition [3] states that the differential dimension of 
the prime differential ideal *p is n — r. Following [20^ Chapter II, Section 12, Theorem 6], 
the transcendence degree of the fraction field of the domain Ai/(Ai n over the ground 
field K equals (n — r)(i + 1) + c for all i sufficiently big, where c is a non-negative integer 
constant which depends only on the ideal This constant c is called the order o/*p and 
denoted by ord(*P). 

The polynomial H<^iT) := (n — r)(T + 1) + ord(*p) is called the Hilbert-Kolchin poly- 
nomial of ^p. The minimum of the indices io where H<;p{i) = trdeg/^Frac(Aj/(Aj n ^i)) for 
all i > io is known as the Hilbert-Kolchin regularity of the ideal *p. 



trdegj^(Frac(5i/Ai„e+i+fc n Bi)) = {n - r){i + 1) + er - /ifc. 



(3) 




i=l j 
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4.1 The Hilbert-Kolchin regularity 



A well-known result from the theory of characteristic sets states that the Hilbert-Kolchin 
regularity of a prime differential ideal is equal to max{ord(C) : C G C} — 1 where C is a 
characteristic set of the differential ideal for an orderly ranking (see \20\ Ch. II, Section 
12, Th. 6 (d)] and [1 Theorem 3.3]). 

The results developed so far enable us to exhibit the following simple upper bound 
for the Hilbert-Kolchin regularity of the ideal *p depending only on the maximal order 
of derivation of the variables involved in the system ([2]), independently of knowing a 
characteristic set: 

Theorem 12 The Hilbert-Kolchin regularity of the ideal *P is bounded by e — 1. 

In particular, for first-order systems of type |^ (in other words, for the case e = 1), 

trdegj^Frac(Ai/(^i n = {n-r){i + l) + ord(«P) for all i € Nq. 

Proof. Since for all i G Nq the fraction field of the domain Ai/[Ai n *P) coincides with 
the residual field of the local ring obtained from Ai after localization at the prime ideal 
Ai n we have that Frac(Ai/(^i n = Frac(Sj/(5i n A)). Then, we have that 
trdeg^Frac(Aj/(j4j n*p)) = trdeg;^Frac(i?j/A n Bi) and so, it is enough to show that, for 
aU i > e - 1, trdeg^Frac(5j/A Bi) + n - r = trdeg^Frac(5i+i/A n Bj+i). 

Fix an index i > e — 1. Due to Theorem [TUl we have that A fi i?j = Aj+i^g+o- H Bi 
and A n = Aj_|_2_e+f7 n -Bj+i. Thus, using Proposition [6] we obtain: 



Hence, the result holds. ■ 

We point out that as a consequence of this theorem we deduce that the order of any 
characteristic set of *P for an orderly ranking is bounded by the maximal order e. 

4.2 The order 

As a consequence of the results of the previous sections we are able to prove the following 
characterization for the order of the ideal 

Proposition 13 Let *p C K{Xi . . . be a minimal prime differential ideal contain- 

ing polynomials /i,...,/,. G K{Xi . . . ,Xn}- Assume that the DAE system defined by 
fi, . . . , fr is quasi-regular at Then, ord(*P) = er — fi„. 

Proof. Let io be the regularity of the Hilbert-Kolchin function of Then, for every 
i > if), the polynomial H<;p{i) agrees with the Hilbert-Kolchin function of ^P, that is 
(n — r){i + 1) + ord(*P) = trdeg;^Frac(^j/(*P n Ai)). Localizing at the prime ideal *P H Ai, 
we have that Frac(^i/(q3 n A^)) = Bi/^i = Bi/{A n Bi) for every i > io, where A is the 
differential ideal generated by /i, . . . , in K{X}. 



Since Bi^ is a Noetherian ring, for h G Nq sufficiently big, we have that the equality 
A n Bi„ = Ah n Bio holds and so, (n - r)(io + 1) + ord(*p) = tTdegK{BiJ{Ah D Bi,)). 



tidegKFTac{Bi+i/{B,+i n A)) 
tTdegKFrac{B,/{Bi n A)) 



(n - r)(i -1- 2) + er - fi^, 
{n — r){i + 1) + er — fi^^. 
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Now, by Proposition m trdeg^(Sio/(A/, n BiJ) = {n- r){io + 1) + er - ^h-io+e-i, which 
imphes that ord(*p) = er — i^ih-ig+e-i for h sufficiently big. Therefore from Theorem [8] we 
have that ord(^P) = er — fi^^. ■ 

5 Jacobi-type bounds 

In [16j and [1^, Jacobi introduces a parameter associated to the orders of derivations in a 
DAE system and conjectures an upper bound for the order of the system in terms of this 
number. 

Definition 14 Let A € Nq^"*, s <m, be an integer matrix. The Jacobi number of A is 

defined to be 



Jacobi asserts that the order of a DAE system (S) with r = n equations is bounded 
above by J{£), where £ := {eij)i<ij<n is the matrix whose entries are the maximal 
derivation orders Cij of the variables Xj appearing in the polynomials fi and — oo whenever 
the variable Xj does not appear in /j. 

We consider also the integer matrix £q which consists in the previous matrix £ setting 
instead of each — cxd. Clearly J{£) < J{£q)- 

In this section, we will show that the sum of differentiation index of a quasi-regular 
DAE system at a prime differential ideal ^ and the order of ^ can be bounded in terms of 
the maximum and minimum orders of the equations and the Jacobi number J{£q). With 
the same techniques we can also recover the Jacobi bound for the order of the system (S) 
(see Theorem 1181 below). 

Our main result of this section is the following: 

Theorem 15 Let (S) be a DAE system defined by polynomials fi, ■ ■ ■ , fr S K{Xi, . . . , X„} 
and let ^ he a minimal prime differential ideal containing fi, . . . , fr such that (S) is quasi- 
regular at ^. Consider the integer matrix £q := [Uj)i<i<r,i<j<n, where eij := ordxj (/j) 
and if the variable Xj does not appear in fi. Then, the '^-differentiation index a of the 
system (S) and the order ord(*p) of the differential ideal *p satisfy 



Before proving the Theorem 1 15t we will show some auxiliary technical results concern- 
ing the Jacobi number of integer matrices. 

Let L be an arbitrary field. Given a matrix A = (oij) G L^^*", we can estimate its rank 
by means of the Jacobi number of an associated binary matrix: we define B(A) ^ L^^"^ 
as the matrix whose entries are 



s 



J{A) :=max| ^ai^(i) | r : {1, . . . , s} ^ {1 



m 



} is an injection 



a + ord(*P) < J{£o) + max{ejj} — minjejj}. 
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Lemma 16 Let A € L^""" with s <m. Then rank(A) < J{B{A)). 

Proof. Set J := J{B{A)). By definition, J = maXr{J2i<h<s ^i^)hT{h)} i where the 
maximum is taken over all injections r : {1, . . . , s} — > {1, . . . , m}. Since B(A)ij = or 1 
according as aij equals zero or not, it follows that for every r, the tuple (aiT-(i)i • • • i o,st{s)) 
has at most J nonzero entries. 

Let k with J < k < s and choose indices 1 < ii < ■ ■ ■ < ik 1^ s and 1 < ji < • • • < 
jk ^ nT-- Consider the k x k square submatrix of A corresponding to the intersection of 
the rows indexed by ii, . . . and the columns indexed by ji, - ■ ■ ,jk- The determinant 
of this submatrix equals Z^,^ sg(!/) ni</i<A; 01^,1/(1^)) where the sum runs over all bijections 
f : {ji,---,jk}- Now, for each u, the vector (ui^^i^i^), . . . ,ai^^(^i^)) has at 

least one zero entry (since k > J and every can be extended to an injection r) and so, 
the determinant vanishes. We conclude that rank(^) < J. ■ 

Our second auxiliary result relates the Jacobi number of an integer matrix to the Jacobi 
number of an associated binary matrix of a particular kind. 

Fix € N. For every < a < A;, let „ G Nq^'^ be the lower triangular matrix defined 

by 

_Jl ifk-a<i-j<k-l 
I U otherwise, 

i.e. Tk^a has the a lower diagonals with all their entries equal to 1 and all the remaining 
entries of the matrix are zero: 



a < 



•■■ 






.. \ 











1 


... 






1 ■•• 








1 • ■ ■ 


'•• 

1 1 





..■ 1 



Given a matrix A = {aij) G Nq^'" and a positive integer k > maxjj{ajj}, we define a 
ks X km block binary matrix Tk{A) as follows: 



^'^k,aii • • • -^fc,ai„ 



Tk{A) 



\Tk,asi 



Tk. 



Lemma 17 Let A e Nq^"^ with s < m. Then, for every k > maxj j{cLij}, we have 
J{A) = J{Tk{A)). 



13 



Proof. First, note that it suffices to prove the identity for square matrices: in any case, 
one can make A square by adding null rows, which does not change its Jacobi number or 
the Jacobi number of its expansion Tfc(A). So, we assume A G N™^™. 

By definition, J{A) = ^8,-^TeSm{J2i<i<m^iT{i)}- Without loss of generality, we may 
assume that the maximum is attained when r is the identity. We consider a permutation 
r € Skm which chooses a /cm-tuple of entries of r/c(^) consisting of the an I's of the 
upper nonzero diagonal of each block T^^au in the main diagonal of Tk{A) along with 
some zero elements: for instance, the permutation defined by T((i — l)k + h) = (i — l)k + 
rk{h + Ojj — 1) + 1 for every 1 < i < m, 1 < h < m, where denotes remainder in 

r(0 



the division by k. Then, XlflTi ^fc(^)/f(0 = S2=i "^m — It follows that J{A) < 




max^e5,„ Ei=i Tk{A)if^i) = J{Tk{A)). 

In order to prove the other inequality we will use the following result known as the 
Konig-Egervary theorem (see [10], |23j): 

m m 

J{A) = minj ^ Ai + ^ </> j I >^i,(pj e No, Xi + (t)j > aij VI < i, j < m|. (4) 
i=i j=i 

Consider nonnegative integer vectors A = (Ai,...,Am) and (j) = [4>i^ ■ ■ ■ ■•(km) such that 
Ai + (t)j > aij for every 1 <i,j <m and J{A) = Aj + YlJLi 't>j- 

We define nonnegative integer vectors A = (Ai, . . . , A^), <I> = (<I>i, . . . , € Nq™ 
with Aj, G Ng as follows: for every 1 < i < m and 1 < j < m, 

) for ft, = 1, . . . , A; — Aj 
I 1 for h = k — Ai + 1, . . . , k 

Let us show that A,$ satisfy + ^ji > Tk{A)(^i^i)f,+h,{j^i)k+i, or equivalently, 

A,ft + «5j7 > (5) 

for every I < i, j < m and 1 < h,l < k. 

Recall that {Tj,ai)hi is nonzero only for those ft, / with k — a^j <h — l<k — 1 and, in 
these cases, it equals 1. 

If ft > fc — Aj + 1 or Z < we have Aj/j + <I>j/ > 1, since at least one of the terms equals 
1 and so, inequality ([5]) holds in this case. Now, if ft < A; — Aj and I > — 1, both Aih and 
^ji are zero, but h — I <k — {Xi + — 1 < k — a^j — 1, which implies that iTk^aij)hi = 0. 

Finally, since 

m k m k m m 

i=l h=l j=l h=l 1=1 j=l 

applying identity to the matrix Tk{A), we conclude that J{Tk{A)) < J{A). m 

We are now ready to prove the main result of this section. 
Proof. (Proof of Theorem] 1^ ) We will show that 

er — J{£o) < fJ-k 1^ — ord(*p) V /c > e — minjejj}, (6) 
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where e = max{ejj}. 

These inequahties imply that the sequence {fJ-k)k>e-mm{eij} increases at most (er — 
ord(*P)) — (er — J{£o)) = J{£o) — ord(*p) times and hence the complete sequence {nk)keNo 
increases at most e — mm{eij} + J (Sq) — ord(*P) times. By Theorem [S] this number bounds 
the differentiation index a and Theorem 1151 follows. 

Then, it suffices to prove the inequality 

First, note that by Theorem[8]and Propositionll31 the inequality /u^ < /.to- = er— ord(*p) 
holds for every fc S Nq, which proves the second inequality. 

Now, fix A; > e — min{ejj}. By Defintion[5l fik = kr — vank{^k,e-i) ■ In order to simplify 
notation, we will write 5^ := ^k,e-i- Let us consider the matrix ^fc which is obtained by 
permutation of rows and columns of ^k so that: 



3fc 





-11 


9/f 


-1] 






-11 




ffe-ii 
-1] 






-fc-lj 




ffe-ii 




dfl'' 


-1] 


dfl" 


-1] 








dxi^^^ 


ffe-ii 





A 



where for every 1 < i < r and 1 < j < n. 



dfi 



df 



[k-i] 



[e,e+k-l] 



dX. 



df. 



dX 





dfi 



dX 



(e+i) 



df^ 



(fe-i) 



dfi 








dX 



(e+l) 



dX, 



(e + fc-l) 



kxk 



(7) 



Applying Lemma [TBI we have that rank(5/fc 
going to estimate the Jacobi number of i?(3A:)- 

As ordxj (/i) = ^ij, the partial derivative 



dfi 



dX 



rank(5fc) < J{B{2k)). We are now 
is identically zero for all p, q with 



(p) 



q — p > Cij. Thus, the only nonzero entries of the block 



df, 



■fc-ii 



dX 



[e, e+fe-1] 



or, equivalently, all 



the I's of the corresponding block in the binary matrix lie in places {h,l) with 

h — I > e — eij. 

Taking into account that the Jacobi number of a matrix does not decrease by replacing 
some of its zero entries by positive integers, we deduce that the Jacobi number of B{^i.) 
is bounded above by the Jacobi number of the block binary matrix which is obtained by 
setting to 1 each entry of corresponding to all these partial derivatives that could 

be nonzero: 



k—e+eii ^fe,fe— e+ei2 
^fc,fc— e+e2i ^fcjfc— £+£22 



/Tk 
Tk 

\Tk,k—e+tri '^k,k—e+tr2 



Tk,k-e+.,^ 

Tk,k-e+t2n 

Tk^k—e+ernJ 



(8) 
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Note that the above matrix is the binary matrix T{Sk) associated with the integer 
matrix 

fk-e + eii k-e + ei2 ... - e + ei„\ 

k - e + €21 k - e + €22 ■ ■ ■ k - e + €2n 

\k - e + €ri k - e + €r2 ■ . . k - e + e™/ 

Therefore, due to Lemma [T71 J{B{Zk)) < J{T{Sk)) = J{Sk) = ™axT-{^[^-^ A: — e+ejT-(j)} = 
{k - e)r + max^{X;[=i eir(i)} = {k - e)r + J{8o). 

We conclude that rank[Zk) ^ {k — e)r + J[£Q) and so, = kr — ia,'Bk{Zk) > er — J(£'o) 



Remark on the proof Suppose that the variable Xj does not occur in the equation 
fi\ then the matrix d?]) is zero. Hence, the corresponding triangular block T^ ^.e+e-^ in 
([H]) may be taken as the zero matrix and we can replace the integer k — e + €ij by in 
the matrix always preserving the inequality rank(5fe) < J{T{Sk)) = J{Sk) for the new 
matrix S^. Therefore, if we redefine €ij := — oo for those pairs (i, j) such that Xj does not 
appear in /j, we conclude that the inequality 

r 

iwk{Zk) < max{y~^ max{0, k — e + ej^(j)}}. (9) 

holds for any integer k verifying k > e — minjejj}. 

This remark is a key fact we use to obtain the Jacobi bound for the order of the prime 
ideal *P: 

Theorem 18 Let *p C K{Xi . . . , Xn} be a minimal prime differential ideal contain- 
ing polynomials /i,...,/,. E K{Xi . . . ^X^}. Assume that the DAE system defined by 
fi,...,fr is quasi-regular at ^ and verifies the hypothesis of Subsection \ 2.4\ Let £ := 
{^ij)i<i<r,i<j<n, where €ij := ordxj(/'t) and —oo if the variable Xj does not appear in fi. 
Then, ord(«p) < J{£). 

Proof. It suffices to show that the inequality ([6]) holds for some index k after replacing 
the matrix £q by £. Since the inequality //fc < er — ord(*p) is independent of the matrix £ 
it remains to prove that there exists an integer k such that er — J{£) < holds. From 
the definition of fi^ this inequality is equivalent to rank(3fc) < (k — e)r + J{£). 
From ([9]) in the previous remark it is enough to prove that 

r 

max{y~^ max{0, k — e + ejT-(j)}} < {k — e)r + J{£) 

^ 1=1 

holds for some k > e — minjeij}. 

Take any k > er and let r : {1, . . . , r} — > {1, . . . ,n} injective. 

First, suppose that ejT-(j) / — oo for all i = 1, . . . , r. Then max{0, k — e + eir(i)} = 
k — e + €i^(j^-^ for all i. Hence ^1^=1 maxjO, k — e-\- £^^-(4)} = {k — e)r + ^iT{i)^ which is 
clearly bounded by {k — e)r + J{£). 
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Now, suppose that ei()T(io) = ~°o for some index iq. Then we have 



^max{0, k-e + eir{i)} = ^ max{0, k-e + ei^^i)} < k{r - 1), (10) 

because each max{0, A; — e + eir(i)} is bounded by k. The quasi-regularity of the sys- 
tem ensures that J{£) ^ — oo (see for instance [22]); then there exists an injection 
: {1, . . . , r} — > {1, . . . , n} such that ei,^(i) 7^ —00 holds for alH = 1, . . . , r. Then 



{k — e)r < max{0, k — e + ei^(i)} < {k — e)r + J{£) 



(11) 



1=1 



holds. Since k > er, we have k{r — 1) < (k — e)r. Hence, from pop and (jlip we have that 
max{0, fc — e + ej^(j)} < k{r — l) < {k — e)r < {k — e)r + J {£) . The theorem follows. ■ 

The following two classical examples of DAE systems show that the bounds in Theo- 
rems [T2] and [TH] may be attained. 

Example 1. Let us consider a system in Hessenberg form of size n (see, for instance, [H 
Definition 2.5.3]): 

-^1 = fl{Xi, . . . ,Xn) 

(S) := { X, = /,(Xi_i,X„...,X„_i) 2<i<n-l 



(12) 



^ = UXn-i) 

such that (^ gx{-i ) ■ ( axnia ) " " " (^) " (^) ^ assume that the differential 

ideal *p := [Xi - /i, . . . , Xn-i - /n-i, /n] is prime. 
The matrix £q associated with this system is 



/ 1 


. 


. 


\ 





1 . 


. 








. 


. 1 





Vo 


. 


. 


/ 



£0 



and we have the following upper bound according to Theorem [T5l 

a + ord(*P) < J{£o) + maxjejj} — minjey} = n — 1 + 1 — = n. (13) 
In order to compute the differentiation index of this system, consider the matrices 



A :-- 



( 1 


. 


. 


\ 





1 . 


. 








. 


. 1 





V 


. 


. 


/ 
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and 



B. 



a,f3 









a/1"' 


\ 




axf' • ■ ■ 




ax 


(f) 






a/("' 













sxf' ■ ■ ■ 






















• • ■ 














d.fi"'> 







/ 


ding matrices ^^,0 (see Definition H]) 


are 


A 













\ 




A 












-62,1 


B2,2 


A 










for a, /3 > 1. 



\ -^(fc-i),i -fc^(fc-i),2 



(fc-l),3 



It can be shown that the dimension of the corresponding kernels of Z\ g ^'^^ A*fc = ^1 
for A; = 0, . . . , n. Therefore, the differentiation index of the system (S) is cr = n and the 
order of the differential ideal *P is ord(^) = 1-n — /i„ = n — n = 0. Thus, the upper bound 
(fT3|) is attained. 

In order to verify the upper bound stated in Theorem 1181 we consider the matrix 



1 


★ 




-k 










1 








—00 


—00 





1 


★ 




—00 


—00 


—00 


—00 . . 


1 




—00 


—00 


—00 


—00 . . 


. 


1 


—00 


—00 


—00 


—00 . . 


. —00 





— CXD 



\ 

where each -k stands for or —00. 

It is easy to see that the Jacobi number of this matrix is J{£) = 0, which agrees with 
the order of 



\ 



Example 2. Consider the DAE system arising from a variational problem describing 
the motion of a pendulum of length L. If 5 is the gravitational constant and A the force 
in the bar: 



X. 



(2) 



(2) 



- \Xi 
\X2+g 



Xf + xl - L2 



(14) 



Let <p := - XXi,X^''' - XX2 + g, Xf + X'i - L'^] c M{Xi, X2, A}. 



.(2) 
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2 0\ 

The matrix £q associated to the system is I 2 I and the bound given by Theorem 

VO 0/ 

US] is 0- + ord(*p) < J{£q) + max{eij} - min{eij} = 4 + 2- = 6. 

The corresponding matrices 3fc,i) = 1, . . . ,n (see Definition H]) are 



1 









































^ 





1 































































































1 



































n 
u 


n 
u 


n 


n 


1 
i 


n 

u 


n 
u 


n 
u 


n 
u 


n 
u 


n 

u 


n 
u 


n 
u 


n 

u 


u 















































-A 





-^1 











1 





























-A 


-X2 














1 























2Xi 


2X2 









































-3A 





-3X1 


-A 
















1 




















-3A 


-3X2 





-A 


~X2 














1 














QXi 


6X2 





2X1 


2X2 
































-6A(2) 







-4A 





-4X1 


-A 
















1 











-6A(2) 







-4A 


-4X2 





-A 


-X2 














1 





12Xf' 


12Xf' 





8X1 


8X2 





2X1 


2X2 
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The dimension of the corresponding kernels of Z\ i are //^ = A;, for k = 0, . . . , 4, and 
/i5 = 4. Then, the differentiation index of the system (S) is o" = 4. We remark that 
the pendulum system is customarily considered as a system of differentiation index 3, but 
in fact this is the value of the differentiation index of its first-order equivalent system; 
as we observed in [8l Subsection 5.1] differential changes of coordinates may alter the 
differentiation index. 

On the other hand, the order of the associated differential ideal is ord(*p) = er — fi„ = 
2-3-4 = 2 (see Proposition [Tl). 

Therefore, we have that cr + ord(*p) =4 + 2 = 6, which coincides with our stated upper 
bound. 

/ 2 -00 \ 

Finally we have ^- = I — 00 2 I . The Jacobi number of this matrix equals 2 
\ -00/ 

and thus, the Jacobi bound in Theorem 1181 coincides with the order of the ideal. 



6 The membership problem for quasi-regular differential 
ideals 

Roughly speaking, the ideal membership problem in an arbitrary commutative ring A 
consists in deciding if a given element f ^ A belongs to a fixed ideal I C A, and, in 
the affirmative case, representing / as a polynomial linear combination of a given set of 
generators of /. 
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Since the work of G. Hermann [T3], it is well known that the membership problem is 
decidable for polynomial rings in finitely many variables. On the contrary, this is not the 
case in the differential context, where the problem is undecidable for arbitrary ideals (see 
|12| ) and remains still an open question for finitely generated ideals. However, there are 
special classes of differential ideals for which the problem is decidable, in particular the 
class of differential radical ideals ([M], see also [2]). 

Concerning the representation problem, besides the non noetherianity, the differential 
case involves another additional ingredient: the order N of derivation of the given gen- 
erators of / needed to write an element / € / as a polynomial linear combination of the 
generators and their first N total derivatives. The known order bounds seem to be too 
big, even for radical ideals (see for instance [13], where an upper bound in terms of the 
Ackerman function is given). Obviously, once the order is bounded, the problem becomes 
a purely algebraic representation problem instance in a suitable polynomial ring in finitely 
many variables. 

By means of Theorem [10] above, we are able to give efficient order bounds for the 
membership problem in the quasi-regular differential setting which lead to estimations for 
the degrees and number of variables involved in the representation. The quasi-regularity 
condition also ensures that the (algebraic) ideals involved are prime and generated by 
regular sequences, which allows a significant improvement in the degree bounds for the 
representation with respect to more general situations. 

We will use the following membership theorem for polynomial rings: 

Theorem 19 (f^ Theorem 5.1]) Let k he a field and gi,...,gs £ k[Zi, . . . , Z^] be a 
complete intersection of polynomials whose total degrees are bounded by an integer d. Let 
g E k[Zi, . . . , Zn] be another polynomial. Then the following conditions are equivalent: 

1. g belongs to the ideal generated by gi, . . . , gs,' 

2. there exist polynomials ai, . . . , such that g = Y1 ^j9j '^^^ deg(aj(7j) < d^ + deg{g) 
for 1 < j < s. ■ 

With the same notations as above we have the following effective differential member- 
ship: 

Theorem 20 Let (S) be a quasi-regular differential system in the sense of RemarklM and 
Section \K^ defined by polynomials F := fi, . . . , fr € K{X}, and let D be an upper bound 
for the total degrees of fi, . . . , fr. Let f £ K{X} be an arbitrary differential polynomial in 
the differential ideal [F]. Set N := a + max{—l, ord(/) — e}, where a is the differentiation 
index o/(S). Then, a representation 



holds in the ring Aj\f^e where each polynomial gij has total degree bounded by deg(/) -|- 




l<i<r 
0<j<N 
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Proof. The upper bound on the order of derivation of the polynomials fi,. . . ,fr is a 
direct consequence of Theorem [in] applied to i := max{e — 1, ord(/)}. The degree upper 
bound for the polynomials gij follows from Proposition [3] and Theorem [T9l ■ 

Remark 21 From Theorems\T^ and\T^ we deduce that for every i € N>e-i, the equality 
Aj_|_i+j(£-y)„inijj{^.^.} n Bi = A n Bi holds, which provides a completely syntactical upper 
bound for the derivation orders and degrees in the membership problem: it suffices to take 
N := J(fo) — min{ejj} + max{ord(/), e — 1}. 

7 DAE and ODE systems 

7.1 The ^-differentiation index and an explicit ODE system 

In the zero-dimensional case {n = r), the estimation for the Hilbert-Kolchin regularity of 
the ideal ^ allows us to give a result concerning the number of derivatives of the input 
equations required to obtain an explicit ODE system from the system ([2]). We will show 
that this number is at most the ^^-differentiation index of the system: 

Theorem 22 Let (S) be a DAE system as in (0j of differential dimension (or equiv- 
alently, r = n), maximal order bounded by e and ^-differentiation index cr. Let H = 
■ ■ ■ jCs^"^} ^ {^['^"-'^1} be an algebraic transcendence basis of the fraction field of 
-Be-i/A n -Be-i over the field K . Then: 

1. for each i = 1, . . . , s there exists a non-zero separable polynomial Pi with coefficients 
in the base field k, such that ^|^^) G {f[^\ • • • , /i*^') C B^+a,' 

2. set . . . , r/„} := {X} \ {^i, . . . , ^s}. Then, for all i = s + 1, . . . ,n, there exists 
a non-zero separable polynomial Pi with coefficients in the base field K , such that 

P.(H,^f-^))G(/f-^l,...,/['^-^l)ci3e+.-i. 

In particular, for every i = l,...,n there exists a separable non-trivial polynomial 
relation between X^^'^ and H modulo A which can be obtained using at most a derivations 
of the input equations. 

Proof. Since we are in a differential zero-dimensional situation, from the upper bound 
on the regularity of the Hilbert-Kolchin function (Theorem I12|). we have that the set S is 
also an algebraic transcendence basis of the fraction field of the ring B^/A D B^. Then, 
for i = 1, . . . , s, there exists a polynomial Pj in s -|- 1 variables with coefficients in K, such 
that belongs to the ideal A D Be = A^+i n Be (Theorem [10]) . Clearly, this 

polynomial can be chosen separable. 

The second assertion follows similarly, but in this case we use the fact that the family 
{rj^'^ ^\ H} is algebraically dependent when regarded in the fraction field of Be-i/AnBe-i 
over K, for all i = s -|- 1, . . . , n. ■ 
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7.2 Existence and Uniqueness of Solutions 



In the case of (explicit) ODE systems, it is well known and classical that, under certain 
general hypotheses, one can ensure the existence and uniqueness of solutions when certain 
initial conditions are fixed. However in the case of implicit DAE systems the same problem 
has only been considered recently (see for instance [32], [57]). 

This last subsection is devoted to the problem of the existence and uniqueness of 
solutions for an ubiquitous class of zero-dimensional implicit autonomous DAE systems. 
More precisely, throughout this subsection we consider DAE systems 



where /i, . . . are polynomials in the n differential unknowns X := Xi, . . . ,Xn and 
their first derivatives, with coefficients in C. We also assume that the equations generate 
a prime zero-dimensional differential ideal £}:= [/i, ...,/„] C C{X} such that the system 
is quasi-regular at and satisfies the hypothesis of Section 12.41 We denote by a the 
differentiation index of the system (6). 

Under these assumptions (cf. Remark [2] and Proposition [3]) , the ideal Qa+i generated 
by the first a total derivatives of the polynomials fi is a prime ideal defining an algebraic 
variety W contained in C"^°""''^\ 

Set Vi C C^" and Vq C C" for the projection of W to the coordinates {x, x) and 
X, respectively. So, Vi and Vq are defined by the contraction of the prolonged prime 
ideal l2o-+i to the corresponding polynomial rings C[X, X] and C[X] respectively. Let us 
remark that from Theorem \T0\ these varieties are also defined by the contractions of the 
ideal O. C C{X} to the corresponding rings. 

Roughly speaking our existence and uniqueness theorem will state that for any point 
{xq,pq) in a suitable Zariski open dense U Vi the system admits a unique solution 
if : {—e,e) — > C" such that ip{0) = xq and ip{0) = pq. As customarily (see for instance [571 
Section 6]), the determination of the open subset U needs additional properties concerning 
the smoothness of both the initial point (xo,Po) and the related projection tt : Vi ^ Vq, 
(x, x) 1-^ x. 

In order to formalize these ideas we recall the notion of unramifiedness (see for instance 
[251 §6, page 100]). 

Definition 23 Let S be an R-algehra, p C S a prime ideal and q C R its contraction to 
R. The ideal p is unramified if qSp = pSp holds and the field extension of the residual 
fields k{p)\k{q) is separable algebraic. 

If F : V ^ W is a morphism of complex algebraic affine varieties and p V , we say 
that F is unramified at p if the maximal ideal associated to p is an unramified prime 
ideal ofC[W] for the structure of C[V]- algebra induced by F (observe that in this case the 
separability condition is automatically fulfilled). 

We recall that the unramified points form an open Zariski dense subset of the source 
space (see for instance [251 Corollary 6.10]). 



f fiiX,X) 







(6) := 



(15) 



I fn{X,X) 
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Let us reinterpret the condition of unramifiedness for our projection vr : V^i ^ Vq at a 
point (xo,Po) £ Vi: if 9Jt C C[Vi] denotes the maximal ideal associated to {xo,po), then 
the contracted ideal dJl n C[Vo] agrees with the maximal ideal C C[Vo] associated to 
the point 7r(xo,Po) = a^o ^ ^o- So, the unramifiedness says that 9^1 must generate the 
maximal (hence regular) ideal of the local ring C[Vi]fnt- Equivalently, we have the equality 
9Jt = I{Vi) +91 in the local ring C[X, X]fxn (here I{Vi) denotes the ideal of the variety Vi). 

Therefore, if H is any system of generators of I{Vi), the polynomials H{xo,X) gen- 
erate the maximal ideal of po in the polynomial ring C[X] localized at the maximal ideal 
corresponding to pq. Hence, in terms of the Jacobian matrix DH{xq,po), the unramified- 
ness of vr at the point (xo,Po) is equivalent to the fact that the submatrix of DH{xq,pq) 
corresponding to those derivatives with respect to the variables X has full column rank n. 

Let us remark that under this last form the unramifiedness condition becomes very 
similar to the conditions induced by the hypotheses of ^^non- singularity" and ''^reducible 
TT-manifold" given in Rabier & Rheinboldt's paper (see |37^ Theorem 5.2 and Definition 
6.1]) required in order to prove their existence and uniqueness result \37\ Theorem 6.1]. 
On the other hand, our requirement is less restrictive than the birationality condition 
assumed in [351 Theorem 25]. 

Now we are able to show the main result of this subsection: 

Theorem 24 Let (6) be a DAE system as in fl^). Let 0. = [/i, . . . , fn] C C{X} be the 
associated differential ideal and satisfies the hypothesis of Section \2.4\ Assume that O is 
prime and (©) is quasi-regular. Let Vq C C" and V\ C C^" be the irreducible varieties 
defined by the ideals Q R C[X] and £2 H C[X, X] respectively and let tt : Vi ^ Vq be the 
projection to the first n coordinates. 

Then, for every regular point (xo,Po) £ where vr is unramified, there exist e > 0, a 
relative open neighborhood lA dVi of {xq^pq) and a unique analytic function ip : (— e, e) — > 
C" such that {ip, p) : (— e, e) — > U, ip{0) = xq and for every i = 1, . . . , n, fi{ip{t), p>{t)) = 
for t G (-e,e). 

Proof. Let G C C[X] and H C C[X,X] be systems of generators of O n C[X] and 
n C[X, X] respectively. Without loss of generality, we may assume that G C H. 

From the unramifiedness condition the submatrix D^^H consisting of last n columns 
of the Jacobian matrix DH has full column rank n at (xo,Po)- In addition, since {xq,pq) 
is a regular point of Vi, the rank of the matrix DH{xq,pq) equals 2n — d, which is the 
codimension of Vi. Then, without loss of generality (after a permutation of the variables 
X), we may assume that the submatrix of DH{xo,po) consisting of its last 2n — d columns 
has full column rank. 

Let us write x := Xi, . . . , X^ and y := X^+i, ■ ■ ■ , Xn for the variables X, and xq := 
(xo,yo) G X C"~'^. By the Implicit Function Theorem (IFT), there exist open neigh- 
borhoods Zq C of xq and Zi C C^'^"'^ of (y^jPo) and holomorphic functions {h,^,r]) : 
Zq —5- Zi such that the open subset of Vi given by {H = 0} n (Zq x Zi) can be described 
by the new, explicit system 

y = /i(x) 

X = e(x) . (16) 

y = ??(x) 
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Note that if if = {G, G}, the Jacobian matrix of H at {xo,po) takes the form 



D^G{xo) 


DyG{xo) 





D^G{xo,po) 


DyG{xo,po) 


D^G{xo,po) 



ip{t) = {tp{t), h{ilj{t)) is a solution of the mixed differential system 



DH{xo,po) -- 

Then, we have that rank(DyG(xo)) = n — d and so, by the IFT applied to the system 
G = at the point xq, the system y = /i(x) describes a neighborhood of xq in Vq. 

Let us consider now the system x = ^(x) as an autonomous ODE. By the classical 
Existence Theorem (see for instance \Vb\ Theorem 1, §3, Chapter 8]), there exist e > 
and a unique analytic function ip := (-f/;!, • • • , V'rf) '■ ~* C'^, such that -(/'(O) = xq 

and, for j = 1, . . . ^ifjj = o -0. Thus, the function if : {—e,e) — > C" defined by 

y = /i(x) 

X = e(x) ■ 

It remains to verify that 99 is a solution of the original system {&)■ Since (/i, . . . , /„) C 
£1 n Ai = (H), it is enough to check that if verifies the system (fT6]) since it defines an 
open subset of Vi = {H = 0} containing the point (xo,po)- By construction, one only 
needs to check that ^ fd+j = Vj ° ^oi all j = 1, . . . ,n — d, or, more explicitly, that 

For every (7 € H Aq, we know that g{x, /i(x)) = in a neighborhood of xq. Taking 
the partial derivative with respect to the variable Xi, for i = l,...,d, we have that 
^-(x, /i(x)) + ^"~^ dXa+ ^(■^)) ~ ^ X in a neighborhood of xq. Multiplying 

this identity by ^i(x) and adding over all the indices i, we obtain: 

n n—d d rj, 

j=l j=l ■' i=l 

On the other hand, since the total derivative of g belongs to the ideal 0.0 Ai, we have 
that g{x, /i(x), ,^(x), r/(x)) = for x in a neighborhood of xq, that is, 

d n n — d rj 

^-f (x,Mx)) e.(x) + j;^^(x,Mx)) 7?,(X) = 0. (18) 
From identities (|17p and (|18p we deduce that, in a neighborhood of xq, 

Ea^(''."W)(%w-E||w«.w) = »^ 

In particular, by using this identity for all g £ G, we conclude that, for all x in a neigh- 
borhood of xq, the vector w(x) £ C"""', whose j'-th coordinate is '>]ji^)~Yli=i S^(x) Ci(x) 
for j = 1, . . . , n—d, belongs to the kernel of the submatrix DyG{x, h{x)). Since this matrix 
has full column rank n — d, it follows that w{x) = and thus 

dh- 

^i(^) = Yl 'djt^^^ ^'^^^ j = l,...,re-d. 

i=l * 



24 



Setting x = ijj{t) we conclude that 



The uniqueness of the solution is a straightforward consequence of the classical theorem 
for ODE'S. ■ 

We remark that the previous existence and uniqueness result can be easily extended 
to non autonomous systems in, at least, two well known ways. 

First, suppose that the base field K in the square system (jlSp is the field of rational 
complex functions C{t) instead of C. Let r{t) G C{t) be the maximum common divisor of 
the denominators appearing in (jlSp and consider the new system 



(Si 



' r{t)h{X,X) = 

; , (19) 
_ r{t)fniX,X) = 



which involves polynomials in C[t, X, X]. 

This system can be interpreted as an autonomous one by simply considering t as a 
new unknown variable and adding the equation 

i = 1. 

It is easy to see that this new system is quasi-regular, the involved equations generate 
a differential ideal in C{t,X} which is prime after inverting the polynomial r{t) (in fact 
the corresponding localization of the factor ring can be naturally included in the factor 
ring of the original system which is an integral domain by hypothesis). Also the system 
verifies the hypothesis of Subsection 12.41 as the original one. 

Moreover, the localization in r(t) by the Rabinowitz's Trick can be replaced by intro- 
ducing a new variable S and the (polynomial) equation 

Sr{t) -1 = 0. 

Now, we can applied the previous result to this (n + 2) x (n + 2)-autonomous system 
and extend Theorem 1241 for the non-autonomous case. 



On the other hand, as it is observed in [361 Section 2.2], we may also extend our results 
to non-autonomous systems involving transcendental functions which may be defined as 
solutions of single DAE-equations (for instance, trigonometric or exponential functions): 
one replaces the transcendental function by a new variable and adds the autonomous 
differential equation which defines that function. Since each added equation involves a 
new variable the assumptions about the original system will be fulfilled also by the new 
one. 
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